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I.  INTRODUCTION 


A  set  of  FORTRAN  subroutines  for  contour  plotting  has  been  prepared. 

These  subroutines  accept  arrays  of  X,  Y,  and  Z  coordinates  as  input  and  draw 
contours  of  constant  values  of  Z.  The  present  version  has  been  prepared  for 
use  oh  the  Interdata  system  located  in  the  Systems  Simulation  and  Development 
Directorate  of  the  US  Army  Missile  Command  (MICOM)  at  Redstone  Arsenal, 

Alabama.  The  subroutines  can  be  adapted  to  any  other  computer  installation 
having  output  graphics  capability  by  rewriting  certain  subroutines  which 
provide  the  commands  to  the  output  graphics  devices.  A  later  section  of  this 
documentation  will  outline  those  changes  necessary  to  transfer  to  another 
system. 

The  source  program  for  the  complete  set  of  subroutines  is  contained  in 
Appendix  A.  THe  source  program  is  in  FORTRAN. 

II.  OPERATION  AND  FEATURES 

The  contour  plotter  operates  by  finding  the  intersection  of  a  plane 
parallel  to  the  X-Y  plane  with  a  plane  constructed  so  as  to  pass  through 
three  neighboring  grid  points.  The  theory  of  this  operation  is  presented 
in  the  Theory  Section.  Since  the  intersection  of  two  planes  is  a  straight 
line  segment,  the  contours  are  comprised  of  straight  line  segments.  Clearly, 
the  smaller  the  mesh,  the  more  contour  detail  will  be  revealed.  It  is  not 
necessary  that  the  X-Y  grid  be  uniform  or  even  rectangular. 

In  normal  operations  the  subroutines  will  automatically  carry  out  scaling 
operations  to  set  up  appropriate  scale  factors.  When  this  mode  of  operation 
is  selected,  the  scaling  will  be  such  as  to  include  the  maximum  and  minimum 
values  of  X  and  Y. 

If  desired,  it  is  possible  to  enable  an  option  which  causes  the  plot 
to  "zoom  in"  to  cover  any  desired  rectangular  region.  In  this  case,  the 
plot  scales  are  automatically  adjusted  so  that  the  desired  region  fills  the 
display  space  on  the  output  graphics  device. 

As  previously  noted,  the  grid  does  not  have  to  be  rectangular.  It  should 
be  noted  that  sometimes,  the  use  of  non-rectangular  grids  may  result  in  plots 
which  display  a  skew  distortion. 

A  contour  numbering  option  is  available  which  causes  the  contours  to  be 
numbered  sequentially.  If  the  numbering  option  is  used,  it  is  necessary  to 
specify  a  parameter  which  identifies  the  number  of  points  on  a  contour  to  be 
plotted  between  the  numbering.  If  the  number  of  points  between  numbering 
is  given  as  zero,  the  contour  will  be  numbered  the  first  time  that  a  point 
is  plotted  and  not  thereafter.  Not  all  contours  need  be  numbered. 

The  contour  plotting  package  is  not  limited  in  the  number  of  grid  points 
which  may  be  supplied  or  in  the  number  of  separate  value  of  Z  for  which  contours 
may  be  drawn.  However,  it  must  be  noted  that  large  arrays  of  grid  points  may 
make  the  plotting  process  very  slow.  In  the  present  version,  the  maximum 
number  of  contours  which  may  be  labeled  is  fifty. 
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III.  THEORY  OF  OPERATION 


The  basis  for  operation  of  the  contour  plotter  lies  in  several  elementary 
geometric  concepts.  Specifically: 

•  The  intersection  of  two  planes  is  a  straight  line 

•  A  plane  can  be  constructed  which  passes  through  any 
three  distinct  points  in  space 

•  A  straight  line  segment  is  uniquely  determined  by  its 
end  point 

In  actual  operation,  the  plotter  works  by  scanning  through  the  grid 
in  a  columnar  manner  considering  four  neighboring  grid  points.  A  typical 
set  of  four  neighboring  grid  points  is  shown  in  the  figure  below.  This  set 
of  points  forms  the  complex  A,B,C,D.  In  this  figure,  upper  case  symbols 
represent  points  in  the  three  dimensional  space  X,  Y,  Z  and  the  lower  case 
symbols  represent  the  projection  of  the  corresponding  upper  case  point  onto  the 
X,  Y  plane. 


In  order  to  plot  a  contour,  the  continuous  surface  in  the  region  of  the 
complex  is  approximated  by  the  four  triangular  plane  areas  (or  platelets) 
shown  in  the  figure  as  ABG,  BCG,  CDG,  and  DAG.  Points  e,  f,  and  g  all  lie  on 
the  intersection  of  the  diagonals  ac  and  bd.  Thus,  the  points  E,  F,  and  G  all 
have  the  same  X  and  Y  coordinates  as  the  intersection  point  g.  These  coordinates 
are  easily  computed  from  the  coordinates  of  the  projected  points  a,  b,  c  and  d. 
Thus, 
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<XCVA  -  XAYC><XD  -  V  -  <XDVB  -  XBYD)(XC  -  V 
^  ‘  <YD  '  V  (XC  -  V  -  <YC  -  V  (XD  -  XB> 


<xayc  -  xcya)(yp  -  V  -  (Vd  -  V.X'c  -  V 

(XD  -  XB)(YC  -  V  -  <XC  '  XA)(YD  -  V 


Using  linear  interpolation,  the  Z  coordinates  of  points  E  and  F  are  computed 
as: 


h  •  h  *  <zc 


ZA> 


(xG  -  xA)2  4-  (yg  -  ya)2 

<XC  -  XA>2  +  (YC  -  YA>2 


ZP  •  ZB  +  (Z0 


ZB> 


(xg  -  xb>2  *  (vc  -  v2 
«D  -  V2  +  <YD  -  V2 


The  Z  coordinate  of  G  is  chosen  to  be  the  mean  of  Z  and  Z  ,  i.e., 

E  r 


z 


G 


(ZE  +  V 
2 


Each  of  the  four  platelets  (e.g.,  ABC)  is  examined  to  determine  if  a  plane 
parallel  to  and  a  distance  Z  above  the  XY  plane  will  intersect  with  that 
platelet.  In  order  to  make  this  determination,  the  Z  values  of  the  vertices 
of  the  platelet  are  placed  in  ascending  order.  The  desired  Z  value  is  then 
compared  to  the  largest  and  smallest  values  of  the  vertex  values.  If  the 
desired  value  lies  within  the  range,  the  plane  will  intersect  the  platelet. 
Further,  if  the  desired  contour  value  is  greater  than  the  intermediate  point  Z 
value,  the  contour  plane  will  intersect  the  edges  of  the  platelet  which  connect 
the  vertex  having  the  maximum  Z  value  with  the  vertex  of  minimum  value  and  with 
the  vertex  of  intermediate  value.  If  the  contour  value  is  less  than  the 
intermediate  point  Z  value,  the  contour  plane  will  intersect  the  edges  of  the 
platelet  which  connect  the  vertex  having  the  maximum  Z  value  with  the  vertex 
of  minimum  value  and  with  the  vertex  of  intermediate  value.  If  the  contour 
value  is  less  than  the  intermediate  point  Z  value,  the  contour  plane  will 
intersect  the  edges  of  the  platelet  which  connect  the  vertex  of  minimum  Z 
value  with  the  vertex  of  maximum  value  and  with  the  vertex  of  intermediate 
value. 

The  coordinates  of  the  edge  intersection  points  are  then  determined  by  a 

linear  interpolation.  Let  (X  ,  Y  ,  Z),  (X  ,  Y  ,  Z  ),  and  (X , 

max  max  max  lnt  int  int  min 

^min*  ^min^  t*le  or<*ere<*  sequence  of  vertices  of  a  platelet.  If  Z  > 
then. 
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r 


(z  -  z  ) 

v  _  y  .  _ max 

' 'i  ‘  max  (Z  -  Z  )  ''''int  ''max' 

int  max 


(X,  -  X  ) 


Y  =  Y  + 


(Z  -  Z  ) 
max 


(Y,_  ~  Y  ) 


i  max  (Z  -  Z  )  int  max 
int  max 


(Z  -  Z  ) 

y  —  Y  _i_  max 

A.  “  A  T  /  v  “  / 

j  max  (Z  .  -  Z  )  min  max 

min  max 


(X_,^  -  X  ) 


Y.  =  Y  + 


(Z  "  Z  aJ 

max 


(Y  .  -  Y  ). 


j  max  (Z  .  -  Z  )  min  max 

min  max 


If  Z  <  Z.  then 
int 


X  =  X  .  + 


(Z  -  Zmin> 


-  X_,  J 


'i  min  (Z,  -  Z  .  )  v  int  min" 

int  min 


Y  =  Y  ,  + 


(z  -  z„i„> 


a...  -  t_,j 


i  min  (Z.  -  Z  )  int  min 

int  min 


X,  -  X.  + 


<Z  -  Z,ln> 


(X_  .  -  X  ,_) 


j  min  (Z  -  Z  )  max  min 
J  max  min 


Y  =  Y  .  + 


<z  -  w 


(V _ -  V  J_). 


1  min  (Z  -  Z  .  )  max  min 
J  max  min 


A  line  is  then  drawn  from  the  point  (X^,  Y^)  to  (X^ ,  Y^) , 


This  process  is  carried  out  for  all  four  of  the  platelets  at  the  end  of 
this  cycle.  At  this  point  a  portion  of  the  contour  is  complete  as  shown  by 
the  line  mnop  in  the  figure.  The  entire  process  is  repeated  for  each  set  of 
four  neighboring  points. 


It  should  be  noted  that  a  somewhat  similar  algorithm  could  be  designed 
using  only  two  platelets  to  represent  the  surface.  In  this  case,  a  diagonal 
would  be  drawn  between  two  of  the  diagonal  points.  For  example,  the  platelets 
might  be  ABC  and  ACD.  This  choice  of  platelets  sometimes  results  in  a 
skewing  of  the  contours.  The  four  platelet  scheme  essentially  eliminates  this 
skewing.  To  the  best  of  the  author’s  knowledge,  the  four  platelet  scheme  has 
not  been  previously  described. 
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IV.  CALLING  SEQUENCE 

To  utilize  the  package  as  a  contour  plotter,  the  following  FORTRAN 
statement  must  be  used: 


CALL  SLICE 

<ZV,X,Y,Z,M,N,i 

MM, NN,NZ, MS, ZOOM, XLO, XUP, YLO, YUP, TYPE, TAG, 

NTAG 

,TAGV,NTIMES) . 

The  arguments  in  this 

call  have  the 

following  meanings  and  characteristics. 

Argument 

■TXBl 

Dimension3 

Definition 

ZV 

REAL 

NZ 

The  Z  values  of  the  contour  to  be 
drawn 

X 

REAL 

MMXNN 

The  X  values  of  the  grid  points 

Y 

REAL 

MMXNN 

The  Y  values  of  the  grid  points 

Z 

REAL 

MMXNN 

The  Z  values  of  the  grid  points 

M 

INTEGER 

The  number  of  "columns"  in  the  X-Y 
grid 

N 

INTEGER 

M 

The  number  of  grid  points  in  each 
"column"2 

MM 

INTEGER 

The  maximum  number  of  "columns"  in 
the  grid3 

NN 

INTEGER 

The  maximum  number  of  "rows"  in  the 
grid3 

NZ 

INTEGER 

The  number  of  contours  to  be  drawn 

MS 

INTEGER 

If  the  zoom  option  is  not  selected 
and  if  MS-1,  automatic  scaling  of 
the  grid  will  occur. 

ZOOM 

LOGICAL 

If  ZOOM  is  true,  the  ZOOM  option  is 
selected 

XLO 

REAL 

Minimum  X  value  of  the  plotting 
window 

XUP 

REAL 

Maximum  X  value  of  the  plotting 
window 

YLO 

REAL 

Minimum  Y  value  of  the  plotting 
window 

*If  not  a  simple  variable. 

2Each  "column”  need  not  have  the  same  number  of  points. 

3Variable  dimension  statements  are  utilized  for  several  variables. 
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Argument 

Ty£e 

Dimension1 

Definition 

YUP 

REAL 

Maximum  Y  value  of  the  plotting 
window 

TYPE 

INTEGER 

N2 

Sets  the  type  of  line  to  be  drawn 
for  each  contour.  Permissible 
values  are  fLINEf ,  rDASHf ,  'DOT' 

TAG 

LOGICAL 

If  TAG  is  true,  the  contours  specified 
in  TAGV  will  be  labeled  with  their 
contour  number  the  first  time  a  point 
is  plotted  in  the  window  and  every 
NTIMES  points  thereafter. 

NTAG 

INTEGER 

The  number  of  contours  to  be  tagged 

TAGV 

INTEGER 

NTAG 

A  vector  containing  the  contour 
numbers  of  the  contours  to  be  tagged2 

NTIMES 

INTEGER 

The  number  of  points  plotted  within 
the  plotting  window  between  labeling 
of  the  contours.  If  NTIMES  -  0,  the 
contours  will  be  tagged  at  the  first 
point  plotted  only. 

The  first  use  of  subroutine  SLICE  or  any  of  the  other  subroutines  which 
utilize  output  graphic  pen  motion  must  be  preceeded  by  the  following  two  calls: 

CALL  INIT 

CALL  DEVICE  (NDV) 3 

Successive  calls  to  SLICE  should  be  separated  by 
CALL  PAGE 

which  erases  the  screen.  A  call  to  subroutine  COPY 
CALL  COPY 

causes  a  hard  copy  to  be  produced. 

1 If  not  a  simple  variable 

2Not  all  contours  need  be  tagged  if  the  tag  option  is  selected. 

3The  device  number  NDV  is  dependent  upon  the  output  graphic  device  being 
used:  NDV  *  1  for  HP2648  Terminals 

NDV  *  2  for  TEK4006  and  TEK4051  Terminals 

NDV  *  3  for  TEK4014  Terminals 

NDV  *  4  for  TEK  4662  Flat  Bed  Plotter 


If  this  package  is  to  be  adapted  to  another  computer  system  or  graphics 
package,  it  will  be  necessary  to  rewrite  certain  subroutines.  These  subroutines 
provide  basic  graphics  operations  such  as  pen  motion,  etc.  In  the  paragraphs 
to  follow,  these  subroutines  and  their  functions  are  described. 

SUBROUTINE  UWINDO  (XLO,XUP,YLO,YUP) 

This  subroutine  should  set  the  upper  and  lower  bounds  on  the  X  and  Y 
coordinates  of  the  display  window.  These  values  are  respectively  XUP,  XLO, 

YUP,  YLO.  Note  that  these  values  are  in  virtual  units,  i.e.,  number  being 
plotted,  and  are  not  meant  to  be  physical  distance  units  on  the  display. 

SUBROUTINE  UMOVE  (X,Y) 

This  subroutine  causes  the  pen  to  move  from  its  current  location  to  a 
point  having  virtual  coordinates  X,Y.  This  pen  movement  is  accomplished 
with  the  pen  "up”,  i.e.,  no  line  is  drawn. 

SUBROUTINE  UPEN  (X,Y) 

This  subroutine  causes  a  line  to  be  drawn  from  the  current  pen  position 
to  the  point  having  virtual  coordinates  X,Y.  If  any  portion  of  this  line 
occurs  outside  of  the  display  window,  that  portion  of  the  line  should  be 
visible . 

SUBROUTINE  UPRNT1(N,B) 

This  subroutine  causes  an  integer  number  N  to  be  printed  at  the  current 
pen  location.  If  that  point  is  outside  of  the  display  window,  the  printing 
should  be  invisible.  The  argument  B  is  a  dummy. 

SUBROUTINE  USET(A) 

This  subroutine  is  used  to  define  the  line  type  that  is  to  be  drawn  for 
all  subsequent  calls  of  UPEN  until  changed  by  another  call  of  USET.  If 
A  =  'DASH*,  the  line  is  a  dashed  line.  If  A  «  'LINE1  a  solid  line  is  drawn. 

If  A  -  'DOT1  a  dotted  line  is  drawn. 
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subroutine  sliceizv.xx , yy,zz.m,*#mm. nn.nz, ms, zoom. xlo.xup 

YLO. YUP.l YPF. TAG,NTAG,TAGV,NTIMFS) 

LOGICAL  ZOOM. TAG. nFwTAG 

DIMENSION  XA<*m.Nn),YY<MM,NN).ZZ(MM.NN).N(M),X(3),Y(3),Z(3) 
,XPC3),YR<3)  ,  ZP  (  3  )  , 2 V ( NZ )  ,TX(4)#TYM),TZ(4),NEi*TAG(50) 

INTEGER  TYPE  (n/.)  ,TAGV  (nTAG)  , MTAG(50) 

♦*♦**4 *♦***•*♦♦**♦♦***♦*♦• ♦♦♦♦*****£***#%*♦***♦*♦*♦♦♦••♦•♦♦*•♦*♦♦•♦*♦ 

THIS  SUBROUTINE  DRA*S  A  CONTOllk  VaP  TOR  THE  2  COMPONENTS  * 

CORRESPONDING  TO  A  A  SET  OF  MESH  POINTS.  THE  MESH  DOES  NOT  * 

MAVt  TO  bt  RECTANGULAR.  IP  IT  IS  DESIRED  TO  SEf  UP  APPROPRIATE  * 
SCALE  FACTORS.  vS  SHOULD  BE  1  * 

ZV  IS  ^  ARRAY  CONTAINING  THE  VALDES  OP  Z  FOR  WHICH  CONTOURS  * 

ARE  DFSlKfcU  * 

If  ZOO**. TRUE.  .  The  FIELD  OF  VlFw  WILL  he  zoomed  in  on  the  * 

RECTANGULAR  A  H  fc  a  SPECIFIED  bY  XlO , XUU . YLO , YUP .  IF  .NOT. ZOOM  * 

XLO.XUP# rtO.iUP  *IlL  RETURN  THE  iNFATi  VALUES  COMPUTED  FOR  * 

OSCaLF  * 

KM  6NO  Hn  are  IhP  **  A  X l mum  DIMENSIONS  OF  THE  DOUBLE  SUBSCRIPTED  * 
VAhiAbLES  AS  olvK'J  IN  THE  ORIGINAL  ABSOLUTE  DIMENSION  STATEMENT* 
THF  I N I  FGEk  vtaoH  i  1Y PEI  IS  USED  TO  SET  THE  LINE  TYPE  FOR  EACH* 
CiiMOuR*  Ik  ThF.  Ll.vfc.  TYPES  ARE  NOT  SPECIFIED,  THEY  DEFAULT  TO  * 
SOLID  LlNtS.  THE  DASH  LENGTH  AND  SPACE  LENGTH  ARE  SET  AT  .0366  * 

inches,  rtpes  Available  include  'line* , 'dot, ♦dash* .  * 

IF  •TAG'S. r^uE.  ,  I  HE  CONTOURS  SPECIFIED  I N  T AGV  WILL  BE  TAGGED* 

^iCh  IrtFjR  C'ihliJUK  HUMBER  THE  FIRST  TIME  THE  CONTOUR  APPEARS  * 

i.J  THE  -UNCmH.  TnF  iUj^nt M  OF  CONTOURS  TO  BE  TAGGED  IS  NTAG  * 

THE  «AaI*U*’  OF  CuHVtS  THAT  CAN  BE  USED  WHEN  THE  TAG  * 

UPJlO#  IS  »S  F  l,  E  C  T  F  D  IS  !>(»•  * 

r he  Tagged  cuntuuns  *jll  be  marked  the  first  time  that  a  * 

PJ|.«T  0*1  f  H  A  r  CONTOUR  IS  PlO  i  TED  wlTHIN  THE  WINDOW  AND  EVERY  * 
niHtS  U^iiS  iMpWK AFTER.  * 

¥  FUSION  7. F  7/ 10 /MO  * 

IFC  •  NOT. TAG  t  Go  lu  1S3 
DEL  A*  •  UOS*  u  i  K'U-  ALU  J 

i) Eu  Y  *  #  •» 0  7  B  ♦  A  ms  l  Y  •  •  P •  Y  LO  ) 

!)  J  1  S  1  (si,*/, 

m  1*  A  G  (  I  Y  s  J 

lb  1  u e  w i ag  c  t j*. false. 

D'J  15/  isl.-a'AG 
15/  f-Ewf  AGC  TAG*  (  i  )  )s.TM>A  . 

153  C'JPillNUF 

m J  ISO  151,  •-?- 

1  F(  1  YPFI 1  J  .  *F.  'LINE*  .Ah),  TIPtm.NE.  'DASH1  .  AND  »  TYPE  (  I  )  .NE.  '  DOT'  > 
l  fPfe*( l )s'LlNE' 

lbU  CD- riNUE 

CALL  UP6ET< 'StTOASh1 .12. ) 

I F { . N u l .ZOUM.Af  D.KS.FU.I  )  CALL  CSC A LF ( X X . Y Y , M , N , M M , N N *X LO , X Up # 
YLO. YuP  > 

IF(ZUOM)  CALL  0<lMHi(XLO,XUP.YLU.YUP) 

MF,  AX»M-1 
DO  lOu  1  *  J  f  m  «*  A  a 
fllXA*s,n  (  i  )  -1 

IF('«(XM  I  .I.T.m  1 1  )  I  NMAXsNUti 
Du  1UU  J3i,NMAA 
JM  •nCiT./OO*  >  Tti  51'J 
X  ljsXX  (  I  ,  J  J 
y usy / ( i # j  ) 

XL*XX ( 1 # J ) 

YL=YY(I,JJ 
UO  bUO  K 1  *  1  ,2 


i'1*  b  Mi 

X  i  i  i  ♦  i\  \  -  1 

•  J  =  J  ♦  r  <!  -  1 

AL  =  XX  (h  i  ) 

I  h  l  A  M  '  1  ,  <  *•  I  A  1.1  S  X  X  C  K  I  .  K  J  ) 

I  t  l  i  *  (  '  i  ,  r- .  1  j  .  L  •  .  »  i  J  t  i*st  r  IK  i  #*ti) 

bO"  i  r  t  *  I  l  *  l  •  KJ  I  •  ».  !  .  S  ')  V»»a»  V  ik  I  .i».  j  | 

IM\«  .‘.I.AI,  -K.U,M  .YM^HR.JfL.Gr.yUP) 

.  .A-  |  V 

5 1 1  i.M  <  IsuU,  •] 

ui/)a*ni*i..ii 
i  a  ji  =  uuM,j+i  i 
(  \  (  4  I  S  a  X  i  *  ,0  +  1  ) 
n  ft  » s  *  m  1 1  •.»  i 

'•' j  l  2  J*Y  f  (  t  ♦  »  ••»  ) 
ry  (  *  jsi  V  l  i  +  i  ..»♦  I  ) 
i'¥i*i)sikll,.i*1  » 

'*<  (  l  /  =  i.ain') 

'  7t2)*:SZl  Ul  ,J) 

7*,  l  i  is/,/.  I  i  4  j  ,  J4  1  ) 

i  /  t  *)«?■  /  1 1  #  0  +  J  J 
Ut'sJJlU  1-1X12) 

x  n  =  i  xi 3)*  i  *  i  n 
» « v  s  r  i  c  * )  •  r  i  ( / 1 


:'i>.nt»»si  «2*X3  l  -  f  31  **42 

XC=U  U(i)«nu  >-i’X(  n*lYC3)  >*X<i2-CTX(4>*TY(2)-rxm*TY(4)  ) 
*A Jl ) /UtMitf 

yC»UrXC2)»r<(4l-l,Xl4|4lYl2)  k*Y31*(TX(t  )  *T  Y  (  3  )  -  IX  {  J  )  *TY  (  1)  ) 
*Y4?)/u£*U- 

ZCtsTZd  JfiU(3»*U(1  )  MSUM*(  (  (XC-rxcn  >**2*(YC-TY(1))**2)/ 
(X3l«Xil4|J|4Y3m 

ZC2*m2)*(TZ(  *)-lZ(2)  )*SL*r(  t  (  XC-TX  (  2  )  )**2+(  ¥C-TY(  2)  )**2)/ 
(X42*A<*2+Y42*Y<m  I 
ZCs(£Cl*4C? )/2. 

4( l )=ZC 
X ( l )=XC 
Y  (  l  )  s  Y  C 


00  <*0  0  is(^ 

Z ( ?  )  =  i ZCL  J 
Alb) 
k  (  2  J*  't  ra; 

iFlL.tO.4)  O'  T w*  4.)  I 

iiij)*ua*i ) 

X  f  3  jsl  M  L+  l  J 

Y  A  3  )s?Y(L+l  > 

Gu  TO  4J* 

4‘>J  2,1  4)3)  2U  3 

a(3)»»x(U 

Y  l  3  )s’l  Y(  1  J 

402  GALL  bUtOil  X  ,  k  .  I,  XP  ,  YP.ZP  ) 
uj  iuo  k  =  \,.,z 

IMTYPb(h)  .Eu.  'I)UI  •  )  CAliL  UPSET  (  *  SfcTUASH  •  .  9  2  •  ) 

CALL  1 iTtk(ZVlH )*XPtYP*ZPiTlPE(K)#TAG,NF*TAG(K)#K*NTlNE5»#TAC(K) 
»l>f  LX, DEL Y  , XLU, XUP, YLO. YUP,ZUGM  ) 

1M  I  YPtfC*  )  .(.Cl.  'OUT*  )  CALL  UPSET!  1  5FTDASH  •  ,  1 2  .  ) 

2  oo  cot-  n  woe 

400  CONTINUE. 

100  CONTINUE 
RETURN 
END 

SUBROUTINE  I NlER ( ZV , X , Y , Z  fTYPE , TAG , Nf>TAG , 4 TAG  $ WPTS • MPTS# 
D£LX#US;LY,XLU,XUP,YL0#  YUP-ZOOM) 
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* 


LOGICAL  Tag. NE*T4G,KLIH. ZOOM 
1  L  FEGEK  i\VF. 

♦*♦******♦*♦**♦*#♦♦*♦♦****♦♦♦♦■♦♦♦♦♦****♦*♦**♦********♦**•♦♦***♦♦****♦♦ 

*  THIS  SMtiHUH  r ± AflC  is  CALLED  HY  SLICE  AND  ORArfS  THE  LINE  REPRESENT* 

*  IVG  lHt  l»/rgi*$FCTitjrt  uf  A  HORIZONTAL  PLANE  AND  THE  TRIANGULAR  * 

*  plane  region  * 

*  U-RS1U-.  J.i  7/08/So  * 

♦♦*****♦♦**  ********************** 

M.  i  ^  =  ,r  ALSr,  t 

T*-  (  t  1  .Ort.JtV  .GT.i(  J)  )  GO  To  100 

i  *•  l  h  l  J  I  .  L*  .  L  (  M  )  Go  TO  100 

l.'(  i  YPE.Hr  ,  •oor*  )  (.ALL  JSETCTVPM 

ir  iTKHr.t,  •.m'u  '  i  C A l . I  USM’C  MUSH*  ) 

t  ►  (  £  l  3  )  •  L ."  »  /  ,  (  2  )  t  vio  3  0  2  1  u 

1  ►  (  1 1  2  l  .  Lr.  #  /.  (  l  )  >  viL-  'p:  2  00 

2»'J  X  4  s  (  l  A  l  J  j  ~  Z  v  )*M  1  )  ♦  (  7.  v  •  £  1 1 )  )  *  A  (  3  )  )  /  (  Z  l  J  )  ~Z  (  1  )  I 
Msi  (Zi3i-Zv  )*>  (  l  U{ZV-Z(  l>  )*Y(  3)  )/(Z(3  l-Z(l)  ) 

1 C  l  T  A  o  •  A  .4 o  #  *»  £  »  T  A G  )  GO  fu 
2b 2  CALb  ^JVUH.yp 

1K(ZMG1  .61  ?)  3  GU  10  y«.' 

X5*UZ(2)*ZHn(UtUv2U))*xrS))/(2(2)-7.a)) 
ib=(  U(  2  )-?.i  )*V(l)*<ZvZ(1  )  J*»  (  2)  J/(Z12)*Z(  t  ) ) 

GO  1 u  3 u 

20  X5s(  UU  )-Z<M*x  (  2  )♦(  Z*-ZU)  J*U  3)  )/(Z<3>-Z(2)) 

Y  b*  t  U  l  3  3  «Z  ‘/  J  ♦  y  (  /  )  ♦  (  z  v  -7. 1 2  >)  *  Y  f  3  ) )  /  (  Z  (  3  )  -Z  (  2  J) 

Jo  IEUUlN  l  CALL  CL  IPIX4,Y4,X5,  Y5 , ALU . XUp ,  YLO  ,YUP,KLtP) 

It  (  .NUT. Mil**)  CALL  UPKMX5,  Y5) 
lU‘t  CONTlsJE 
RETURN 

20*  IEUMLE.4U11  GO  10  201 
Gu  TU  203 

210  JF(Z<  .GE.Z(  3  )  )  r, u  TU  iu! 

GO  TO  203 

201  X43ACU 
Y 4  =  Y ( t  ) 

XSsaC2  ) 
i  5  =  Y  (  2  ) 

GU  TO  2 <i 2 

3  ol  M  =  XC2) 

Y4sY( 2  ) 

X5=X(3) 

Yb=Y( 3) 

202  Call  OfrPfc£( X4 , i 4 ) 

IKZOUP)  CALL  CL1PIX4, Y4.X5. Y5.XL0.XUP. YLO, YUP, KLIP) 

X  f  l .GUI • A  Li E )  CALL  UPEN ( Xb • Y  5 1 
RETURN 

250  IF (X4.GK.XL0.AN0.X4.Lb.XUP. AND . Y 4 . GE . Y LU . A 10* Y4 . LE. YUP) 

GO  lo  251 
GO  TO  2  5 1 

251  IfCNPlS.EG.l.OR.HPTS.IfO.NPTSl  go  TO  253 
NPTSs^PTS+l 

GO  TO  252 
2b j  *PVS*2 

IE(LPTS.EO.u  )  TAG*. FALSE. 

CALL  USETI »LlNE‘ 1 
CALL  UH0*E(X4, Y4-DELY) 

X5SX4 

Yo*Y4-UELY 

X7«X4 

Y7«Y4*0ELY 

CALL  CLIPUb,  Y6.X7,  YT.XLO.XUP.YLO.YUP.KLIP) 

IF(.NUr.KLlP)  CALL  UPEN(X7,Y7) 
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call  unovecxa-oelx. Y4j 

XS*X4-DELX 

Y7*Y« 

Yb  =  Y  4 

X7»X4+DELA 

CALL  CLlP(Xb.Yb#X7.Y7,  XLO  .XUP.YLQ.YUP.KLIP) 

IF  ( , NUl,*LlP)  CALL  UPEWCX7.Y7) 

CALL  UPRNTl  (NTAG,  •  1  NTEGF.RNUMBEfc  #  I 

IFCTYPE.t'U.  »OuT'  .OR. TYPE. EU.  ♦OASH*  )  CALL  USETUOASH') 

If (TYPE.EQ. ‘DASH* )  CALL  UPSET ( 1  SET DASH  1  •  1  2* ) 

IFCTYPE.EO. »OOT' )  CALL  UPSET f * SFTOASH 1 . 92. ) 

GO  TO  252 
F.ND 

SUBROUTINE  SO* I  3 ( X X . i Y ,ZZ. X, Y.Z) 

DIMENSION  XXI 3)  .U(3>.ZZ(3)«X(3).Y(3).Z(3) 

♦♦********♦♦♦*♦*♦♦♦ 

4  THIS  SUBROUTINE  IS  CALLED  BY  IMER  AND  SORTS  A  SET  OF  THREE  4 

4  POINTS  INTu  ASCENDING  ORDER  BASED  On  THE  Z  VALUES  * 

»  VERSION  1,0  5/16/80  4 

********* *4+*»4***«***4*****«**«4*****444*****4***44 44 444 4444444444*** 

DU  31  1*1.3 
XU)*XX<X) 

YU  )»YYll) 

31  Z(i)*ZZU> 

DO  2  0  =  1  #  2 
DO  1  1=1.2 

lFCZUJ.LT.ZCI  +  i))  GO  TO  i 
UUmxZU  ) 

DLNXsaU  ) 

Oo*Y*'l  (  I  ) 
zinszn  +  n 
yi  n*i  ( t  +  n 

XU)  =  Xl  !♦!  ) 

2» C  I ♦  1  J*UUc< 

A  (  I  ♦  1  )=i)UMX 
K(Ifl  )*0U*V 

1  CONTINUE 

2  IFlZU).LF#Z(2l.A4D.Z(2).Le.Zm  J  GO  TO  3 

*KlTMb.4)X(j  ),*U),Zin.X(2).Y(2).Z(2).XCJ),Y(3).Z(3) 

4  FuR « At (  /  /  #  *  ERROR  DETECTED  IN  SUBROUTINE  SORT1 //J(3(2X«£15«*Jsj  j 

3  RETuR4 
E.vu 

SUBROUTINE  CSC ALE (X#Y .M.N.MM.NN.XLO* .xup.ylon.yupi 

DIMENSION  X(  MM.txN)  .Y(MM.NN)  ,N(») 

4  THIS  SUBKUUTINE  IS  CALLED  BY  SLICE  AND  SETS  UP  THE  PLOTTING  4 

4  *  I N  D  0  *  IN  VIRTUAL  SPACE  4 

4  VERSION  3,0.  6/24/80  A 

Y*AX*M  •  tfoo 
01)  1  1*1. M 

IM  Y<l.N(  J  n.Gf.YMAX)  YMAX*Y( I.NCIll 

1  CONTINUE 
YmI N*YMAX 
DU  2  1  *  1 . m 

1FCYU.I  i.LT.YMm  YM  IN*Y  ( I  #  t  1 

2  CONTINUE 
XMIN*X( 1 . 1  ) 

AMAX*X(R, 1 j 

CALL  NEATCXRAX.AMlN.XUP.XLOM 
CALL  NEATCYMAX.YMlN.YUP.YLOw) 

CALL  U w IN 00 ( XLOw . XUP . VLOtf.YUP) 

RETURN 


16 


if.  A  f  (  A  -  *  X  ,  > ;  **  J  ^  ,a  1.  P  ,  X I  U>.  > 
»*♦***< 
i  i  srAi.t  a-  i> 


tui* 

SlihWUf;  I  I  Mb 

This  iiUUKUU  r  |Nfc‘  JS  CMLiitf> 

*uk  Pontile  scAbts 

l  •  r> 

i  h  I  A  06  l  A  M  4  *  )  .  L  i’  .  *A  06  (  X  M  J  i  ) 
c;Xs4LM(  A  LI  Mil  i.(ah61mMaJ)J 

2  K=io.*^h:x 

XUPe  l  A  l  i  (  <-AA/MM.  »*•- 
I  fU*.  1  i  XaM  *  I  a 

XI.  0*  SI  A  l.«l<  <  •  j  »/M  l-»  .  )  *rft 

X  LO  A  50  . 

k  r.  i  t  J  w  n 

1  rxsAj'*  ri4ufiUlo(  x*if-«  ?  ’ ) 

GO  TO  ^ 

P  -IO 

6»U4kUM  1*1  flit.  U*l«0i1(  A  1,0.  X UP , YOn • Y UP ) 

L’i-  ft  I  ivU  (  1  M 

£  i  •  *  i*  C  *4  /  h  L  t'  ?  /  *  S  C  ■  *  i  j  t  ,  A  /  t  k  f  i  S  v  ‘  A  |  t  *  Y  7.  P  W 
OjflPU’v  /r  Lr  S/  "  1  O 

X6CALP.S C  w  1  *'>('/)- a  I  .•■«•  1  1  )  J  /  t  xup  -Xi  r»  i 

K  SCALtsi  *  J  1  )  -  •  t-  ».»l  J  )  )  /  (  »  -  y  I 

XZKW*»  j  aLI  1  )  -  A  l.ij 

U  »* 

ke  r  i '  h  < 


I 'VfcS  1  N h.  A  T  I 


P.M.i 

Si»Uk/»UJ  M,t  » 

CAUL  LrfAPH 
r  u  1,1  P  L . 1 1  l  S  , 

:<h  \  \  H 
Pm  I 

O  '  J  rS  | .  i  H  I  I  J  P. 

C  A  LL  »j  U*  I  l  \  , 
(\u»L 

K  ►  I  I  .'  H  ‘ 


>'  M  A  ,  r> 


7  .  1 


.  (  <  , 

J  .  I 


».  O 

S*ip»H'  M  I  1  ■*  c  :Ph-  I  I  (<i,iM 

lM  •  •  O  /.‘I.M/  :  •  *■*  I  #  I  P  (  »  L  M 

111  VM'6I  •»-.  H  >bK  (  *>  1 

r  a oi,  Alpha 
C  A  o l-  oti  r 

fHj'Tr.  I  L 1 1  f  1  0  ti  )  i 

1  On  P»ik*  A  I  {  L  l  ) 

r.r  l  ik  i  ■ 
r  ■  * ' 

I  U'lll  I  »c  t'n*r  i  i  ■'  ,  r'  l 

p  T U K  i. 

p  ii» 

I  I  *  p  I  ■*  > 

1  y  i  <«.e  v.  • 1  »  i  i  .ti.i,  i  i'  rn  (4) 

i  p  (  * .  r  * .  *  t*  u  '  )  C*.  i  i.  * .  4 .1  v  e  i  ^  ) 

I  p  C  •  r  4 •  '  I <  I  >  M  *  ILI.  I.  I  •  l  Y P <  1  ) 

k  r.  t  i-M'  '■ 

p.  * ' 

i»i  h'HU  I  i  .r  1  v.  \  \  (  A  4  ,  Y  •*  ,  A?>,  f  b  ,  XM.#XUP»  Y  1.0  ,  YMP,*blP  ) 
t  <  «L  l  C  ».  I  MJ  ♦' 
h  i ■  1  r  *  •  P  H-'ir.  . 

J  t  i  x  4.<*t  .  A»  I  •  A  #  .  P.  •  A  U  p  •  A  *  I  .XS.GK.XLJ.A^U.Xb.LR.XUP, 

^A.,»  .  /A.'.t,  y  I  •.  A  M’.  H..,r  .  Yur.**M.  .  tb.GP.  Yl.u.  AmU.  Yb.  LE.  YUP) 

«.  b  r  I  Ui<«. 

1MM.LT  .A!,- *. A  U-  1 1« 
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If (X4.GC.XUF.AwO.Xb.GT.XUP)  GO  TO  999 
If t Y 4.L1 . 1 LO. At u. I5.LT. YLO)  GO  TO  999 
If  l  Y4.GT. YoP.AnD.Y5.GT.YUP)  GO  TO  999 
If  (XA.LT.ALU.A^D.Xb.GX.XLUJ  GO  TO  JO 
IF (X4.GT.XoP. AAD.Xb.LT.XUP)  GO  TO  20 
IF ( A  S  *  LT. ALu . AwD. A4.GT.XL0 )  Gu  TO  50 
If (A5.Ul.A0F.AND.X4.LT.XUPi  GO  JO  40 
If  ( Y4.LT.<m.AM>.Y5.GT.  YLO)  Gu  TO  50 
If (Y4.br.JUP. AND.T5.LT. XUPi  GO  TO  60 
1F( Y5.Lr.yt0.AN0.Y4.GT.YL0)  GU  TO  70 
If (YS.GT.U.P. AN0.Y4.LT. YUP)  GO  TO  BO 
GO  'll)  200 
99*  KLlfs.  I’PUt. 

kfcTt'KN 
10  XTSALU 

YT*Y44l Y5-Y4  )• ( XLO-X4 ) / ( X5-X4  I 
X  4  *  A  T 

Y  4  =  YT 

GO  TO  100 
20  XTsXUP 

YT=Y4+(  Y5-Y  *  )•  (  <4*ALP)/  (  A4-XS  I 
X  4s  A  T 

Y  4SY  X 

GU  TO  loo 

io  xjsxm 

risibtl i H-YS )♦ (XLU-A5J/ ( X4-A5) 

XSsXT 

i  t>sY  T 
GO  TO  U<) 

40  X'isXOF 

YTs  V  b  ♦  ( i 4-r*> j*( X 5-XUP ) / ( X 5-X 4 ) 

xbsxr 

Y  ^  s  Y  » 

Gu  IU  |o,.. 

5n  YTsyLU 

XTs>  4+  {  X  b»  *'»)♦(  i  L»J*t  4  )/(  Yb-Y  4  1 
A  +  =  >  T 

GO  TO  100 
60  YTsyUP 

X  I‘SA44(  Ab-X4)4(  Y4-IUP)/« Y4-Yh> 

X  4  =  A  T 
y  4=y  i 
GO  TO  100 
7  0  Y  T  =  Y  LO 

XT  =  Ab*(X4-*S)*n  LU-Y  b 
a5  =  a1 
YbsiT 
GO  TO  lOo 
BO  YTsYUP 

XT=AbF(X4-X5)K lb»YUP)/Clb*Y4* 

X5  =  A  r 
i  b  s  Y  T 

100  JF(  X^.GE.AtU.A\D.X4.t,f  .  AUP.A  ,  l> .  Y  4  .  GE  .  Y  I.C 

4AO0.Y4.tt.Y0P.  AND.AS.UF.XfJ.  A^.1  .Ab.i.F.Xi 

*  Ant).  YS.G&*  YLD.AND.Y5.LK.YuP  I  GO  XU  ?0n 
GO  TO  1 

200  CALL  U*0Vfc(X4,Y4) 

RFTOHN 

Mu 

•  Bfcuti 
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